Photon wave mechanics and position eigenvectors 

Margaret HawtorQ 

One and two photon wave functions are derived by projecting the quantum state vector onto 
simultaneous eigenvectors of the number operator and a recently constructed photon position oper- 
ator [Phys. Rev A 59, 954 (1999)] that couples spin and orbital angular momentum. While only 
the Landau- Peierls wave function defines a positive definite photon density, a similarity transfor- 
mation to a biorthogonal field-potential pair of positive frequency solutions of Maxwell's equations 
preserves eigenvalues and expectation values. We show that this real space description of photons 
is compatible with all of the usual rules of quantum mechanics and provides a framework for under- 
standing the relationships amongst difi^erent forms of the photon wave function in the literature. It 
also gives a quantum picture of the optical angular momentum of beams that applies to both one 
photon and coherent states. According to the rules of qunatum mechanics, this wave function gives 
the probability to count a photon at any position in space. 



I. INTRODUCTION 

The current interest in entanglement and its applica- 
tion to quantum information has rekindled the contro- 
versy surrounding the photon wave function P, [2, S 0) S 
Q . It is still unclear what form a real space photon wave 
function should take, or if one exists. In the standard 
formulation of quantum mechanics, the coordinate space 
wave function is the projection of the state vector onto 
an orthonormal basis of eigenvectors of a Hermitian po- 
sition operator. It has been claimed since the early days 
of quantum mechanics that there is no position opera- 
tor for the photon that allows the introduction of a wave 
function in this way. Contrary to these claims, we have 
recently constructed a photon position operator whose 
transverse eigenvectors form a real space basis. Here we 
will use this basis to obtain a photon wave function that 
is compatible with the usual rules of quantum mechanics. 
We will show that this clarifies a number of previously 
unresolved issues regarding the real space description of 
one photon and multiphoton states. 

In 1933 Pauli fi\ stated that the nonexistence of a den- 
sity for the photon corresponds to the fact that the po- 
sition of a photon cannot be associated with an opera- 
tor in the usual sense. Based on definitions of center of 
mass, Pryce found the k-space photon position operator 
rp = iV — ik/2fc + kxS/fc where Sj are the 3x3 spin 
1 matrices, k is a unit wave vector, and Vj = d/dkj 
Q . This operator does not have commuting components 
which suggests that three spatial coordinates cannot si- 
multaneously have a definite value. In 1949 Newton and 
Wigner sought rotationally invariant localized states and 
the corresponding position operators. They were success- 
ful for massive particles and zero mass particles with spin 
and 1/2, but found for photons "no localized states in 
the above sense exist" @. This result is widely quoted 
as a proof of the nonexistence of a photon position oper- 
ator. It has been proved that there is no photon position 
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operator with commuting components that transforms as 
a vector [IQ] . 

Our position operator [11| has commuting components 
but is not rotationally invariant and does not transform 
as a vector [l2|, and thus it is consistent with the pre- 
vious work. Description of a localized state requires a 
sum over all k, and a localized photon can have definite 
spin in the k-direction, that is it can have definite he- 
licity, but it cannot have definite spin along any fixed 
axis. It is the total angular momentum (AM) that can 
have a definite value along some specified direction in 
space [T^. The position eigenvectors are not spherically 
symmetric, instead they have a vortex structure as is ob- 
served for twisted light [l3|- Compared to the Newton 
Wigner position operators for which transformation of 
a particle's spin and position are separable, the photon 
position operator must incorporate an additional unitary 
transformation that reorients this vortex. 

Valid position eigenvectors cannot violate the 
Hegerfeldt [l^ and Paley- Wiener theorems based 
on Fourier transform theory. Hegerfeldt proved that 
a positive frequency wave function can be exactly 
localized at only one instant in time and interpreted 
this to imply a violation of causality. Bialynicki- 
Birula [13] noted that the Paley- Wiener theorem limits 
g (x) = Jjj°° dkh (fc) exp {—ikx) of the form exp {—Ax^) 
to 7 < 1. He then applied this to separate outgoing 
and incoming exponentially localized spherical pulses in 
three dimensions. However, their sum is not subject to 
the exponential localization limit, as can be seen from 
the form of the /c-integral. Position eigenvectors require 
a sum over all wave vectors, and thus must be a sum 
of outgoing and incoming waves that interfere to give 
exact localization at a single instant in time, consistent 
with the Hegerfeldt theorem. 

Maxwell's equations (MEs) are analogous to the 
Dirac equation when written in terms of the Riemann- 
Silberstein (RS) field vector, proportional to E±icB 
where c is the speed of light in vacuum, E is the elec- 
tric field, and B is the magnetic induction. This sug- 
gests that the photon is an elementary particle like any 
other, and that MEs provide a first quantized description 
of the photon. Use of the positive frequency RS vector 
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as a photon wave function in vacuum and in a medium 
has been thoroughly studied [H, [H, [H, [23|. If a field 
^(1/2) with quantum electrodynamic weighting, fc^/^, is 
used as wave function, a metric factor fc~^ is required in 
the scalar product. The real space squared norm then 
goes as / d^r J d^r'^<-^/'^> (r) • (r') / |r - rf and 

thus its integrand cannot be interpreted as a local num- 
ber density |l1S,j21J. Since the photon has no mass, it has 
been suggested that there is no photon number density, 
only energy density [2^. Photon number density based 
on the Landau-Peierls (LP) wave function 'S'^"^ (without 
the factor fc^/^) was investigated as early as 1930 ^22l. |23| . 
Its absolute value squared is positive definite but it has 
the disadvantage that its relationship to electric current 
density and the electromagnetic fields is nonlocal in real 
space [Ilillll,!!!. 

Returning to field-like ^(^/^^ functions, we will show 
here that it is possible to define a biorthonormal basis 
that gives a local density by combining the eigenvectors 
of an operator with those of its adjoint. This formalism 
has recently been applied to pseudo-Hermitian Hamilto- 
nians that possess real spectra [25] . Such a basis provides 
an interesting alternative to explicit inclusion of a metric 
operator when working with electromagnetic fields. The 
density (r).'S'(-i/2) (j.) is local, but it not positive 

definite since it is not an absolute value squared. Only 
the LP wave function defines a positive definite photon 
density, equal to (r)| . However, we will show that 

the biorthogonal field-potential pair gives the same re- 
sults in most calculations. 

In the present paper one and two photon wave func- 
tions and photon density will be obtained by projection 
onto a basis of position eigenvectors. In Section II the 
photon position operator will be reviewed and the scalar 
product and Hermiticity will be discussed. In Section III 
the orthonormal and biorthonormal eigenkets of the posi- 
tion operator will be obtained in the Heisenberg picture 
(HP). We will then derive photon wave functions from 
quantum electrodynamics (QED) in Section IV by pro- 
jecting the state vector onto simultaneous eigenvectors 
of the photon position, helicity, and number operators. 
We will discuss MEs, photon wave mechanics, and an- 
gular momentum and beams in Sections V, VI and VII 
respectively and then conclude. 



II. POSITION OPERATOR 

We start with a discussion of the photon position op- 
erator. The procedure used in [ll| was to construct an 
operator with transverse eigenvectors of definite helicity, 
cr = ±1. In k-space, it is reasonable to expect that the 
transverse function 

4% (k) = (cfc)" e^^^ exp (-zk . r) /W (1) 

describes a photon located at position r, where Uk = kc 
in vacuum and the parameter x will be discussed later in 



this section. Subscripts denote eigenvalues and Cartesian 
components of the vectors xp and e. Cartesian compo- 
nents are used where it is necessary to avoid confusing 
vector notation. The parameter a allows for both LP 
and field based wave functions. The position eigenvec- 
tors are electric and/or magnetic fields if a — 1/2, the 
vector potential if a — —1/2, or LP wave functions if 
a = 0. This is consistent with the QED based interpre- 
tation that a mode with frequency Uk has energy hcok so 
that the square of the fields gives energy density while 
the wave function gives number density. The spherical 
polar definite helicity unit vectors are 

4"l= (^ + *^0)/V2 (2) 

where 9 and (p are unit vectors in the increasing 9 and 
(j) directions. Periodic boundary conditions in a finite 
volume are used here to simplify the notation, and the 
limit as V oo can be taken to calculate derivatives 
and perform sums. If the wave function ([T]) is a position 
eigenvector it should satisfy the eigenvector equation 

?^"V^l,(k)-r^^l,(k) (3) 

where is the k-space representation of the position 
operator and its eigenvalues, r, can be interpreted as pho- 
ton position. 

The operator arrived at in [11] using the condition ([3]) 

is 

f -HS-k^ cot 6l/yfc (4) 

where 

r^p'> ^ilV -ilak/k + kxS/k, (5) 

is a generalization of the Pryce operator, / is a 3 x 3 unit 
matrix, {Si)^^. = —ietjk, and the component of spin par- 
allel to k, Su = k-S, extracts the helicity a. The opera- 
tor f is essentially the usual k-space position operator, 
«V, with terms added to compensate for differentiation 
of the unit vectors and fc" by V. The term kxS/fc gives a 
transverse vector, while S'k</> cotd/k rotates 6 and 4> back 
to their original orientations, and —ilakjk corrects for 
differentiation of fc". It was proved in [ll| that f*^") has 
commuting components and satisfies the other expected 
commutation relations. 

The photon's position coordinates must be real, and 
this normally implies that the position operator must be 
Hermitian. In the LP case the k-space inner-product is 

^*(0)|§(0)^^^^(0)*(k) $(0) (k) 
k. j 

where (k) and 4'^"' (k) are any two state vectors. 
It can be proved by converting the sum to an inte- 
gral and integrating by parts that /*(o)|f(°)*(°A = 
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^r(o)ijf(o)|^(o)^ which impUes that r(°) is Hermitian. 
The case a = 1/2 with inner-product 



k, j 



was considered in [21} and Integration by parts in 
this case requires differentiation of which again gives 
^*(i/2)|p(i/2)^(i/2)^ ^ ^^(1/2)^(1/2) 1^(1/2)^ ^ proving 

that r(i/2) is Hermitian based on the inner-product (O- 
This leads to the nonlocal real space density discussed in 
the Introduction. Alternatively the inner-product can be 
written as 



^(l/2)|^(-l/2)^ ^Y^^i^m* (k) $(.-1/2) (1^) 
k, j 



by defining ^(-1/2) — ^(^/^^/fc, thus avoiding explicit 
inclusion of the factor and the consequent nonlocal 
real space density. The expectation value of the posi- 
tion operator then satisfies /'S'(i/2) |f (-i/2)^(-i/2) 



venience and the choice x — —nicf) gives [ic 



r(i/2)i$r(i/2)|^(-i/2)^ l^^Q j^pp^y ^l^jg ^j^g localized 

state = ^(") = ■01^?)^ this proves that the eigen- 

value r' is still real. However, the position operators 
r^i/^) and r(~i/^^ = r^i/^^t are not self-adjoint. Opera- 
tors such are these are referred to a pseudo-Hermitian by 
Mostafazadeh [2^. Use of pseudo-Hermitian operators 
and a biorthonormal basis is discussed in more detail in 
the next section. 

In fl^ the position operator was generalized to allow 
for rotation about k through the Euler angle x (^i to 
give the most general transverse basis, 
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It was found that the position operator can be written as 

= D {k^iVk-") (8) 

where D = exp (— iS'kX) exp (— iS's^) exp {—iS2d)- Start- 
ing from a wave vector parallel to z and transverse unit 
vectors x and y, D rotates k from the z-axis to an ori- 
entation described by the angles 9 and 0, while at the 
same time rotating the transverse vectors first to 6 and 
(p and then about k through x- For example, when r^^/^^ 
acts on a transverse field parallel to 4> it rotates it to 
y and divides it by y/uJk, then operates on it with the 
usual k-space position operator iV. It then reverses the 
process by multiplying it by ^^uJk and rotates it back to 
its original transverse orientation. This allows r^^/^) to 
extract the position of the photon from the phase of the 
coefficient of the transverse unit vector. 

The quantum numbers {r, a} index the basis states 
for a given x 0)- The z-axis can be selected for con- 



X - ?y 
2V2 
^ x + iy 
2V2 



(cos 9 — a)e 



i(m(T+l) 



V2 



(cos 9 -\- a) e 



■i(mcr — 1) 



(9) 



For example, x — ~0 (™ — 1) rotates 6 and (p back 
to the X and y axes to give unit vectors that approach 
(x + iay) /\/2 in the 6* — > limit. This is useful in the de- 
scription of paraxial beams since the unit vectors describe 
spin alone. Their coefficients then describe all of the or- 
bital angular momentum so that a factor exp {ilz4>) im- 
plies a z-component of orbital angular momentum equal 
to hlz- 

The spin and orbital AM of a photon are not separable 
beyond the paraxial approximation. For unit vectors of 
the form ([9]) the z-component of total angular momentum 
and photon position operators satisfy [13i] 
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This is just the usual commutation relation satisfied by a 
vector operator and an angular momentum component. 
Here it implies that photon position transforms as a vec- 
tor under rotations about the axis of symmetry of the 
localized states. A photon on the z-axis satisfies the un- 
certainty relation AJ^Arfe ^ 0. Unit vectors of the form 
([9]) contribute a definite z-component of the total AM, 
consistent with {sz,lz\ equal to {— 1, toct -I- 1}, {0,TOcr} 
or {l,ma— 1} with — ma, that is total AM has a 
definite value, but spin and orbital AM do not. 



III. POSITION EIGENVECTORS 

Here we will obtain the eigenvectors of the position 
operators discussed in Section II. The LP form of the po- 
sition operator, r^^\ is self adjoint, has real eigenvalues, 
and defines an orthonormal basis as is usual in quantum 
mechanics. To obtain QED-like fields as eigenvectors, 
the choice a = 1/2 is required. In this section we will use 
the mathematical properties of pseudo-Hermitian opera- 
tors to obtained a completeness relation for the field-like 
photon wave function and investigate how it is related to 
the LP wave function. The operators will be obtained 
in the HP picture, so time dependence as determined by 
the Hamiltonian must also be considered. 

We will start with an examination of the expectation 
values to motivate the use of the biorthonormal formal- 
ism. Any Hermitian operator o satisfies the eigenvector 
equation o |/„) = o„ and its eigenvalues, o„, are real. 
To transform from LP position eigenvectors to fields, 
multiplication by y^uJk is required. Assume that rj is an 
operator with positive square root p = y/rj which will 
equal s/iUk in the present application. We can write 



ifn \o\ f,n) = (/n |P (p ^Op) p ^ \ /,„) = U, 
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where O = p ^dp is a similarity transformation, |^/;^ 



pVlV'n)- The eigenvec- 



p-'\U and |0„) = ((/„|p)t 
tor equation becomes O \Tp„) = On IV'n) ^-^id the eigenval- 
ues and inner-products are preserved. If p is a unitary 
operator, that is if p~^ — p\ then O^^ — O is Hermi- 
tian and |0„) — \4'n) ■ The !?/'„) and |0„) eigenvectors 
are the same, and the usual quantum mechanical for- 
malism is obtained. On the other hand, if p is a Her- 
mitian operator satisfying p^ = p then ^ and 
there are two distinct sets of eigenvectors. We can deal 
with this is one of two ways: (1) The metric operator 
?/ = p^ can be introduced to give the new inner-product 
{<t>n\'n'^4'm) ^"^d work only with the basis. (2) We 
can use the eigenvectors of O and the eigenvectors of 
~ rjOrj'^ which are |i/'„) and |0„) respectively. Since 
(/n|/m) = {fn\pp~^ \ frn) transforms to (0„|V'mL the 
eigenvectors \i^„-^) and are biorthonormal [20|. If 
there is degeneracy, a biorthonormal basis can be ob- 
tained by defining a complete set of commuting operators 
(CSCO). 

The properties of pseudo-Hermitian operators and 
biorthonormal bases have recently been investigated by 
Mostafazadeh and can be summarized as [1^ 

o |v„) = o„ ivj , 6t ^ o: |0„) , 

with the completeness relation 

El^")('^nl= 5] (V„ 1 = 1, 

n n 

where 77 is a metric operator and 1 is the unit operator. 
If p = exists. 



pOp- 



(11) 



is self-adjoint and the eigenvectors 0„ are real. Expec- 
tation values are preserved by the similarity transforma- 
tion, 77. 

To apply this formalism to the photon we take rj = ujk 
and work in k-space. Then o = f'"' is self-adjoint and 
the operators = ft^/^) and O — ff^^/^) have the 

biorthonormal eigenvectors ip'^rl^^ 0^) ^'^^ '^ilj^^'' 0^) 
given by Eq. ([T]) that go as y/uJk and l/^/uJk respec- 
tively as required by QED for the electromagnetic fields 
and the vector potential. The position operators and 
their eigenvectors satisfy 



i'i'i'' (k) = c^i^'i^n 

the biorthonormality condition 



(12) 
(13) 
(14) 



E 



(^W.#rl,)='^'(r-r')<5.,.', (15) 



and the completeness relation 



E 



= 1. 



(16) 



Here is the 3-dimensional Dirac (5-function and we can 
interchange a and —a. The field and the LP operators, 
o, are related as 



at 



1/2^ -1/2 

i^i. OUJi. 



(17) 



consistent with ([8|). This transforms the LP position 
operator f^"^ to r(i/2', introducing an addition term 
—il'k/2k. The momentum and angular momentum op- 
erators Kk and h (— k x iV + S) are unaffected by the 
similarity transformation (I17|) . In the angular momen- 
tum case this is because k x k = 0. 

Time dependence is determined by the Hamiltonian 
H + Hq with 



k,o- 



(18) 



where the zero point terms Hq — ^ huJk/2 which are 
unaffected by the photon state will be omitted here. The 
operator Ok.o- annihilates a photon with wave vector k 
and helicity a and satisfies the commutation relations 

flk.cr, ^k' — Sa,a'S]i.k' ■ The Operators and their eigen- 
kets are time dependent in the HP ^27.] . Using the unitary 
time evolution operator 



U{t) = exp{-iHt), 
the HP position operator becomes 

t^hI^U^ (i) (i) =?^"' 



(19) 



(20) 



) with 



'/'r"^ ) given by Eq. 



with eigenvectors W (t) 

^ in k-space. The coefficient of t in the last term of 
(PO)) is the photon group velocity. 

We can define 1-photon HP annihilation and creation 
operators for a photon with helicity a at position r and 
time t as 



^S., it) - 



k,CTj"k,CT 



V 



(21) 
(22) 



For a — 1/2, Eq. (PT|) implies that the biorthonormal 
pairs are related through 



^(-1/2) 

^(1/2) (i) 



(23) 



analogous to the relationship between the vector poten- 
tial and the electric field in the Coulomb gauge. The 
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1-photon position eigenkets normalized according to ((15)) 
are 



- V'S^i) |0) 



(24) 



where |0) is the electromagnetic vacuum state. The 
projection of ([M)) onto the momentum-helicity basis, 
{|k, cr)}, gives back Eq. ([T]) in the Schrodinger picture. 
The free space operators for a photon with helicity a 
satisfy the r-space dynamical equation 



CTCV X [t) . 



(25) 



The annihilation and creation operators satisfy the equal 
time commutation relations 



E 



(26) 



The Hcrmitian operator describing the density of photons 
with helicity cr, obtained by averaging over the a and —a 
forms, is 

(r, t) = [t) ■ ^['f (t) + H.c. (27) 

The total number operator is 

j d\n^^^ (r, t) = Y, aL.«k,a. (28) 

k,a- 

An alternative linear polarization basis can be ob- 
tained if we define operators that annihilate a photon 
state with polarization in the 6 and directions as 



it) 

4>':' it) 



^Mi it) - <-i it) 



/V2, (29) 



respectively. While the direction of these eigenvectors 
depends on k, they do not rotate in space and time, and 
in that sense they are linearly polarized. The inverse 
transformation is 



i^':! it) 



In free space 





it) 


at 






it) 


dt 



^':\t)+^^^':\t) 



cV X 0J. [t) , 



/V2. (30) 



(31) 



= -cNx^["\t), 



If a = these are the operators introduced by Cook [2J] , 
while if a = 1/2 their dynamics is ME- like. 

The localized definite helicity basis states are eigen- 
vectors of a CSCO, so it is the helicity basis that will 
be used here. Linearly polarized fields can be found by 
taking the sum and difference as in (HH). 



IV. WAVE FUNCTION 

In this section we will obtain one and two photon wave 
functions and photon density by projection onto the basis 
of position eigenvectors found in Section III. This density 
is a 2-point correlation function that is based on the LP 
or biorthonormal basis, rather than electric fields alone 
as in Glauber photodetection theory [28l |. 

The QED state vector describing a pure state in which 
the number of photons and their wave vectors are uncer- 
tain can be expanded as 

I*) = co|0)+^CkX,JO> (32) 

k,o- 

+ ^ E VMc,a;k',<T'Ck,a;k',<T'a^aL,^' |0> + •■ 

where cq = (0|5') , Ck,o- = (0 |ak,o^| , Ck,<T;k',o-' = 
Ck',<T';k,<T = (0 lak.trflk'.fT'l *) > and A/'k^<T;k',fT' = 1 + 
&'kM'^<j.<y' ■ Division by 2! corrects for identical states ob- 
tained when the {k, cr} subscripts are permuted while 
^/H /2 normalizes doubly occupied states. A more gen- 
eral state requires a formulation in terms of density ma- 
trices that will not be attempted here. 

The 1-photon real space wave function in the helicity 
basis, equal to the projection of this state vector onto 



eigenvectors of r^p as {jl'r'',aj 

(r,i)=^ck..e|^^(a;fe)"^:^ 



I* 



ik-r— 'icjfci 



(33) 



where we have used Eqs. (|22j), ([24]) and ([32]). The ex- 
pansion coefficients depend on the choice of basis, for 
example when X — > X + the coefficients Ck.o- — 
Ck,<j exp(— icrAx) analogous to gauge changes of the vec- 
tor potential describing a magnetic monopole in real 
space In any basis the inner-product = 

X^k <7 kk.crl = |ci| where |ci| is the net probability for 
1-photon in state j^*). The free space 1-photon dynam- 
ical equations mirror the operator Eqs. (123)) and (I25p . 
They are 



(34) 



To obtain the 2-photon wave function we can project 
I^P) onto the 2-photon real space basis 



V'r,<T,, it),1pr 



givmg 



(r,r';t,t') 



^r"i., it) V'r"i'.,' it') 



(35) 
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Use of Eq 

r(") 



and 



(t)Vv,l',,' (i')«L"aL,."' |0> then gi 
1 
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k,k" 



to evaluate 



21V 



k,a";k'.cr' 
XCk,CT;k',(T' (^fc'^fc')" 



„(x) ix) ik-r— Mfct ik'-r'— i(^fc/t' 

k,CTj k'jCr'j' 



The bases obtained here provide a real space descrip- 
tion of the multiphoton state that " encodes the maximum 
total knowledge describing the system" as discussed in 
Ref. 3 . The electric field wave function used in 0, Q 
/•3gpr RS vectors in [l], [ll| by themselves do not provide a 
iaasis, and this is the root of the criticism of [2] made 
m The 2-photon wave function (|36p is symmetric in 
agreement with [J, Q • 



_i_„(x) (x) ik-r'-icjfct' ik'-r— icjj./t 
"l"'^k',CT' ,j^k,CT,i''^ ^ 



which becomes a 2-photon wave function if we set t' = t. 
A separate symmetrization step is not required since its 
symmetric form is a direct consequence of the commu- 
tation relations satisfied by the photon annihilation and 
creation operators. 

To obtain an n-photon basis the creation operator can 
be applied to the vacuum n times with each occurrence 
having different parameters r, cr, and j. The state vector 
can then be projected onto this basis to give the n-photon 
term. The result is the symmetric n-photon real space 
function 



{m} 



where 



■t t' 



is a short hand for 



n ut^ I*) (37) 



and m represents the m*'' set of variables, quantum num- 
bers and components {r[™l , t^™! , (jl™] , jl™' }. Generally 
the n-photon states provide more information than can 
be measured. Instead the real space helicity a photon 
density, equal to the expectation value of the number 
density operator, ((27|) . can be defined as 



(38) 
* ) + c.c. 



4")(r,t) - (*|H,(r,i)|*) 



The 0-photon contribution to n is 0, while the 1-photon 
contribution is 

4") (r, t) = i*^")* (r, t) ■ *^-") (r, t) + c.c. (39) 
For the 2-photon state substitution of ([^S)) gives 



nl") (r. 



x*i:<x"L-' (r,r';i,i)-f c.c. 



implying that unobserved photons are summed over. A 
similar argument can be applied to each n-photon term. 
Photons are noninteracting particles and the existence of 
a photon density is consistent with Feynman's conclu- 
sion the photon probability density can be interpreted as 
particle density [29l . [30| . 



V. MAXWELL'S EQUATIONS 

In this section we will show that MEs can be obtained 
from QED in two distinct ways. The first is the conven- 
tional approach of calculating the expectation value of 
operators with all modes populated as coherent states. 
The fields obtained in this way are real and they cannot 
be interpreted as wave functions. The second approach is 
to project the state vector onto the position eigenvectors 
obtained when a field operator acts on the vacuum state 
to give fields proportional to the 1-photon wave function 
components in Section IV. 

If the multipolar Hamiltonian is used, the displace- 
ment is purely photonic, while the vector potential will 
include photon and matter contributions [31[. The vec- 
tor potential operator is a sum over positive and negative 
frequencies, photon and matter parts, and both hclicities. 
We can define 



A(r,t) 
aW (r,t) 

(r,t) 



AW (r,t)+A(-)(r,t), 
A(+) (r,t)+A(+) (r,i) 



(40) 



aW (r,t)+AL7(r,i) 



(+) 



where A'~) = A^+^t and the subscripts m and p denote 
matter and photon parts respectively. The electric field 
and magnetic induction arc then given by 



E = -dA/dt~V(l), 
B = V X A. 



(41) 



In the presence matter of with polarization operator V 
and magnetization M the displacement and magnetic 
field operators are 



D = eoE + V, 
H = B/^o - M, 



(42) 



where SI units are used, cq is the permittivity and fiQ the 
permeability of vacuum, and c = 1/^eo/^o- 

The momentum conjugate to the vector potential is 
— where D i is the transverse part of the displace- 
ment operator |20l . |3]| . These operators satisfy canoni- 

cal commutation relations. Since rpj. ^ and i/Jj. ^ satisfy 
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(l26l) we can choose 



(43) 



This is equivalent to the usual QED expansion of A and 
D and thus is consistent with the operator MEs 



^ dB 
V-B = 0, VxE = - — , 

ot 

V-D = p, VxH=j+^, 



(44) 



where p and j are the free charge and current den- 



sities. In free space Dj^^/-\/eo — 



(1/2) 




where the RS operator is 



/ \/2p.Q as defined in [li 



Coherent states are the most classical, and they can 
be used to establish a connection between QED and the 
real Maxwell fields. Following Cohen- Tannoudji et. al. 
[sH ] the complex Fourier transforms of the classical field 
vectors, 



Vk (t) = j dVV (r, t) 
and the normal variables, 



exp (— ik • r) 



7k W 



£^ it) - ck X 6k (t) 



can be defined. For a coherent state with the pho- 
ton occupancy of mode {k, a} described by the complex 

parameter 7^0-: the average photon number is nk.o- — 
I 1 2 

and the probability amplitude for n-photons is 

exp ^— I 7j(. I ^ ^/\/nI. The quasi-classical coher- 
ent state is a Gaussian wave packet that oscillates with- 
out deformation and with relative number uncertainty 
Ank.cr/fik.cr = 1/ |7kCT|- the limit of infinite pho- 
ton number the electric and magnetic fields oscillate in a 
well defined way as do the solutions to the classical MEs. 
Thus 



-(+) 



coh 



(r,t) = ({7k,.}|A^(+^(r,i)|{7k,.}) (45) 



E 



^k-r— f 



, 2eoUJk 

k,o- 



(46) 



and the fields derived from it behave classically in the 
large photon number limit. 

It is also possible to derive one photon positive fre- 
quency MEs from QED. Since it is still widely believed 



that there is no position basis for the photon, this re- 
sult is new. We can define the I-particle states |Vr.cr) = 

Vy 1.9,0) with V(-) = V(+)t and V(+) = for 

any field operator V such that 



V(+) (r,<) = (g,0|V(+)|vI/ 



(47) 



This can be viewed as the projection of the photon- 
matter state vector state onto the n = 1 term of number- 
position-helicity basis. In the ground state |0) both the 
EM field and any matter present are in their lowest en- 
ergy configurations. The operator V^~' creates 1-particle 
that can be a photon or a material excitation. Since the 
space and time dependence originates entirely in the field 
operators, these functions satisfy ME dynamics. The 1- 
photon MEs are, using 1^^ . 



V • ^ 0, V X E(^ 



(48) 



6>D(+) 



Projection of the state vector onto a basis of 1-photon 
position eigenvectors results in intrinsically positive fre- 
quency electric and magnetic fields defined by (ITf|) that 
satisfy MEs. They can be manipulated to give any of the 
commonly used forms of MEs. 

A wave equation can be obtained from (j48p in the usual 
way to give 



1 9^E(+) 



V X V X E(+) 



9 fdV^^^ 



The terms on the right hand side are the polarization, 
magnetic and external contributions to the time deriva- 
tive of the current density [31]. If there is no magne- 
tization or external current and the polarization is lin- 
ear and isotropic we can write V = eox {k) E which 
can be combined with the S^E^+^/Si^ term. Writing 
e (fc) = eo [1 + X [k)] the angular frequency in (fTB|) is 
LOk — fcci/l -|- X {k). Analogous to the creation of an ex- 
citation of the electromagnetic field (a photon) by D*-^-*, 
the polarization operator T'f") creates a matter excita- 
tion. Energy can be transferred between matter and the 
electromagnetic fields, so the matter and EM modes are 
coupled. Self-consistent solution of the matter-photon 
dynamical equations gives the polariton frequencies ujk 
that determine time dependence. 

The energy, linear momentum and angular momentum 
of the free electromagnetic field are conserved. Their den- 
sities and associated currents satisfy continuity equations 
of the form dp/dt -I- V • j = 0. This can be verified using 
MEs, and the steps in this derivation are still valid if we 
replace the products of classical real fields with Hermi- 
tian linear combinations of products of operators. For 
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example, the current describing the flow of energy den- 
sity ^D(-) • D(+) /2eo + B(-) • B(+) is times the 
linear momentum density 



P(r,i) 



D(-) X B(+) -B(-) X D(+) 



*). (50) 



Together with their associated current densities the com- 
ponents of P also satisfy continuity equations which im- 
plies that / (PrP (r, t) is a constant of the motion. If j^*) 
is a 1-photon state (^-l D^") x B(+) |*) = (^-l D^"' |0) x 
(0|B(+) I*) so that 



P(r,i) = -[D(-)(r,i)xB(+) (r,<) 



+ c.c. 



(51) 



with the fields derived using (gl]), ([H]), and (gT]). For 
a coherent state, the quasi-classical expectation value 
({7k,o-} I 6 x B |{7k,^}) 7^ T>coh X Beoh for small |7k,<,| . 
However ([50)1 can be evaluated exactly using ak,o- |7k o-) 
to give 



1 



P (r, t)^- [T>[-1 (r, t) X B^:^ (r, t) + c.c.J (52) 

with A^J^"* given by (|45p . In cither case the angular 
momentum density is 



J (r,i) = r X P(r,t) . 



(53) 



We are now in a position to compare the classical and 
quantum fields and densities. Eq. (j46p describes real 
fields that are the expectation values for coherent quan- 
tum states. Expectation values do not describe 1-photon 
states since in this case the expectation values of the field 
operators are zero. Instead, it is projection onto a basis 
of position eigenvectors that gives 1-photon positive fre- 
quency fields, proportional to components of the wave 
function. For 1-photon and coherent states momentum 
density can be written as a cross product of fields as in 
(l?T|) and ([5^ . Eq. ([50)1 can be used to interpolate be- 
tween these two extreme cases. 

The density D^^") x B'^+^ can be rewritten as [3l| 



D(-) X B(+) = 



Its first term, equal to 

3 , 3 

2 



D(-) X (v X A(+) 

^d(-VaW-(d(-).v)a 



(+) 



(-1/2) 



is the integrand in the expectation value of the real 
space momentum operator —ifV\7. The last term, 
(D^~) • V) A^+\ also contributes to the flow of energy 



density and has important consequences. It is responsi- 
ble for the spin term in the AM ([55]) . This can be seen 
by writing 



^rx [d(-) -Vj A(+) = D(-)xA(+) 

- (d(-) • v) (r X a(+) 
where a x b = — i (a • S) b gives 

3 



D(-)xA(+) = iyvI/f/^)*;.SvI/(-V2) 
2 ^ ^ ^ 



Since V.D*^ ^ = \ the last term contributes 
J cPrv X pA*-^-* to / (Pr3 (r, t) after integration by parts 
which is zero in the absence of free charge. 



VI. 



PHOTON WAVE MECHANICS 



In this section we will discuss first quantized photon 
quantum mechanics. For dcfiniteness we will refer to the 
Barut-Marlin rules for Schrodinger and Dirac particles 
stated in [s^ as: (a) A basis for the space of wave func- 
tions, which describe all the possible states of a particle, 
is defined by a wave equation, (b) A inner-product is 
defined in the space of the wave functions, (c) Expres- 
sions for the probability density and probability current 
are found. They should form a 4- vector whose diver- 
gence vanishes. The expression for the probability den- 
sity should be positive definite, (d) Operators which cor- 
respond to measurements are defined, in particular, mo- 
mentum and position operators, (e) The eigenfunction 
of the operators, normalized to 1 (in the case of discrete 
spectrum) or a (5-function (in the case of a continuous 
function), are found, (f) The position operator, defined 
in (d), and the inner-product, defined in (b), uniquely 
determine an expression for the probability density. The 
theory is consistent only if this uniquely determined ex- 
pression is identical with the one defined in (c) to satisfy 
a continuity equation. This is a consistency test. 

In brief, these rules apply to the r-space wave mechan- 
ics of a single free photon in free space in the following 
sense: (a) Solutions to p4|) . 



(54) 



(r, t) /dt = CTcV X (r, t) 



include positive and negative frequencies. The negative 
frequency solution can be eliminated on physical grounds 
[Tgl . [29I , thus cutting the Hilbert space in half as is done 
for solutions to the Dirac equation [s^l- (b) The inner- 
product of the wave functions describing states and 



(55) 
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exists and is invariant under similarity transformations 
between a = 1/2 and a = 0. (c) The real number and 
current densities obtained by averaging the a and —a 
densities 



2 / J <^ 



c.c. 



(56) 



satisfy the continuity equation 



5n(") (r, t) 



dt 



V-j(") (r,<)-0. 



(57) 



in the linear polarization basis. If a = 1/2, ([60)) is of the 
form considered by Bialynicki-Birula and Sipe [H, 
while if a = ([6T|) is the form used by Inagaki 
However, (pD|) and (pTjl themselves imply that either the 
helicity or the linear polarization basis can be used in 
combination with field-like a = 1/2 wave functions or 
LP a = wave functions. The operator on the right 
hand sides of ([50)1 and ([M]) is the real space 1-photon 
Hamiltonian. 



The density ieoE • A/h has appeared before in the 
classical context and in applications to beams. Cohen- 
Tannoudji et. al. [sij transform the classical electromag- 
netic angular momentum as 



This can be verified using the wave equation. The 

density n(°) = Y.a ^ positive definite, while J = eo J d^rrx {E x B) 



(62) 



i(i/2) j(i/2)j jg ^ 4-vector that can be written as the 
contraction of second rank EM field tensors with 4- 
potentials. (d) The momentum operator is Hk and the 
position operator is given by Eq. ([4]). (e) The eigenvec- 
tors of these operators are 5-function normalized accord- 
ing to l|15p. (f) The position operator and inner-product 

give the density i ^Vr^jl*) (■0r^"^l*) + cc. Some of 
these points will now be discussed in more detail. 

Both positive and negative frequency solutions of the 
wave equation are mathematically allowed. The classi- 
cal solutions are real, and real waves do not satisfy a 
continuity equation or allow a probability interpretation 
[S^ l . It has been argued by Inagaki for LP wave functions 
that the negative frequency solutions with momentum in 
opposite direction to the wave propagation should be dis- 
carded from the physical photon state [29] . A similar case 
is made by Bialynicki-Birula for elimination of the nega- 
tive frequency fields in field- like wave functions [l^, ■ 

As with MEs the photon wave equations can be written 
in a number of equivalent ways, and this will be consid- 
ered next to allow comparison with the existing photon 
wave function literature. The six component wave func- 
tion 



(a) 
hel 



(a) 
1 

(a) 



in the helicity basis and 



tin 



(58) 



(59) 



in the linear polarization basis can be defined. The 
Schrodinger equation is then, using ([51]) and ([?T|) with 
V X a (S • V) a, 



. 9 / *i' 



(a) 



'dt \ qf^^l 



-iS-V 



iS-V 



(a) 



(a) 



in the helicity basis and 



S-V 

-S-V 



(60) 



(61) 



eo / d\ 



^ Si (r X V) + E X A 



i=l 



by requiring that the fields go to zero sufficiently quickly 
at infinity. Although this looks like an expectation value, 
the fields are classical. In a discussion of optical beams, 
van Enk and Nienhuis [13] separate monochromatic fields 
into their positive and negative frequency parts using 



V = [v(+) exp {-iujt) + V(~^ exp {iojt) 



and obtain for total field linear momentum and AM 



V = -i d 



(+) 



.1=1 

■ 3 



^L»,^+^*(-rxiV + S)yl 



(+) 



(63) 



.(64) 



Here we have assumed the absence of matter in writ- 
ing D — cqE, substituted Af^^^ — jwD(+\ and changed 
the notation a bit for consistency with the present work. 
These are classical expressions, but terms at frequency 
2a; do not contribute to the total momentum and an- 
gular momentum, P and J [35] . They look like the ex- 
pectations values of the linear and angular momentum 
operators that would be obtained using the biorthonor- 
mal wave function pair VeoT^A^^^j^ and 



photon 

The number operator iD*^^' • A'^+^ /2h + H.c. was shown 
previously to be the zeroth component of a four-vector 
obtained by contraction of the second rank EM field ten- 
sor with the four-potential ((/>, A) [36]. This dcmonstates 
that the biorthonormal basis is of value for comparison 
with the existing literature. 

It was noted in Section HI that the biorthonormal 
inner-product is equivalent to the use of a metric opera- 
tor. Using ([^5]) and H — hkc in k-space and substituting 
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H for id/dt the inner- product (ISK)) can be written as 



as in [2l|,l38l. 

The number density is the expectation value of the 
number density operator (|27)) as discussed in Section IV. 
The a = ±1/2 wave function pair gives a real local 1- 
photon density n^^/^\ but this density is not positive 
definite. This can be seen from the following example: If 
l^*) is a 1-photon state that includes only wave vectors 
ki and k2 = ki + Ak both with helicity a where Cki.a = 
Ck2,(T = l/\/2 then 

X cos [Ak • r - (fci - k2)ct)]}/V. 

The cosine term can exceed the spatially uniform time 
independent term due to the -s/fc factors, leading to neg- 
ative values. If ^2 ~ /ci, n^^^^^ is approximately equal to 
the positive definite density, n(°\ however only the LP 
wave function satisfies the positive definite requirement 
exactly. 

It thus appears that LP wave functions are essential 
to a probability interpretation. Field-like wave functions 
can be obtained from the LP wave function by a similar- 
ity transformation, and thus are equivalent to it for the 
calculation of expectation values. The operators given 
by Eq. P^)) in the a = case are identical to the opera- 
tors examined by Cook. The equations that they satisfy 
differ from those for D and B only in that their relation- 
ship to charge and current sources is nonlocal. The LP 
number density has been criticized 0, [20|, [HI, but its 
scalar analog, obtained by taking Fourier transforms of 
the Schmidt modes, has recently been applied to sponta- 
neous emission of a photon byan atom and spontaneous 
parametric down-conversion [s^l . For narrowband su- 
perpositions of plane wave states the distinction between 
the LP and field-like form of the wave function has no 
observable consequences @ . 

The operator (j21|) creates basis states that lead 
to the orthogonal transverse 1-photon wave function 



in the helicity basis. The 

wave function components vl/(-i/2) are proportional to 
the vector potential, while vP^^/^) is related to EM fields. 
Contraction of the second rank field tensor F'^'^ = 
d^Ai" - df^A" with the 4-potential as Fi"'*A^ gives a 4- 
vector [3^. Thus [n 
density is its zeroth component 

In the linear polarization basis the density operators 



(1/2) j(i/2)j jg a^ 4- vector and photon 



n(") (r,i) 



-(Q)t -(-q) -(Q)t -(-q) 

■0r ■''Pr + 4>r 4>r + H.C. 



(65) 



^(a)t ^i-a) ^(-a)t ^ (a) 



with xpj. and 4>^ given by (|29|) . Mandel and Wolf noted 
the convenience of a photon number operator, equal to 

-(0)t -(0) 

rp^ ■ tp^ in the present notation, to the theory of pho- 
ton counting for an arbitrary quantum state '331. Cook 
sought detector independent photon density and current 
operators that satisfy a continuity equation. His opera- 
tors are just ([55)1 if we take a = 0. Inagaki reformulated 
Cook's theory in terms of conventional quantum mechan- 
ics [29]. These authors discuss the restrictions imposed 
by photon nonlocalizability, but the existence of a ba- 
sis of position eigenvectors makes this unnecessary here. 
Our operators describe microscopic densities, and there 
is no restriction based on wave length. 

The Lorentz transformation properties of the a = Q 
photon annihilation operators in the linear polarization 
basis were also considered by Cook 41]. He concluded 
that their continuity equation is covariant in the sense 
that is is related to the field vectors in the same way in 
all reference frames. The Hamiltonian, momentum, an- 
gular momentum, and Lorentz transformation operators 
must conform to the Poincare algebra. Since the position 
operator generates a change in particle momentum, the 
boost operator is closely related to the position opera- 
tor. For a free photon in k-space the Lorentz operator 
corresponding to the a = 1/2 case is [s^ 

K(i/2) = /c(iV) + k X S. 
where K^-^/^) ^ K(i/2)t. Using ^ this gives 

K(") = fc-l/2K(l/2)^l/2 

for the LP boost operator which incorporates the sim- 
ilarity transformation. In k-space this is simple, but in 
r-space it is non-local as discussed by Cook 41]. 

It is stated, in [42] for example, that "the non- 
Hermitian formulation is in most cases a mere change 
of metric of a well posed Hermitian problem. Nonethe- 
less, .. , it has been successfully argued that the non- 
Hermitian formalism is often more natural and simplifies 
calculations." These comments apply here. The choice of 
a does not affect expectation values, the inner-product, 
and the existence of a wave equation and a continuity 
equation. Only the number and current densities them- 
selves are affected. The field and LP bases can be viewed 
as alternative descriptions of the photon state. For most 
purposes fields are more closely related to the physics, 
but the LP basis is needed if the band width is large and 
photon number density is required. 

According to the general rules of quantum mechan- 
ics, for a l-photon state the probability that a photon 
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with helicity a will be found at position r at time t is 

'^a (r, €) . More generally the photon number density 
is the expectation value of the number density operator, 



,(0) 



' (r,t) « nS^/^^ (r,f) given by Glauber [28| de- 

fined an ideal photodector to be of negligible size with a 
frequency-independent photoabsorption probability. An 
ideal photon counting detector also has a quantum effi- 
ciency of 77 = 1 , that is any photon reaching the detector 
is counted. A detector with all of these characteristics 
measures photon position. Consider a 1-photon pulse 
travelling in the z-direction that is normally incident on 
a detector of thickness A2; and area AA. The probability 
that a photon is present in this detector, and hence that 
it is counted, is n'a^ (r,^) AAAz. In Glauber theory the 



E(-)(r,f)-E(+) (r,i) 



count rate is duQ/dt c 

where {duQ / dt) Az /c is the probability the photon is 
counted during the time that it takes to traverse the 



detector. Since 



,(1/2) 



(-) 



a: 



(+) 



Act'*'V^ + c.c where 



-iE^^ /lo for most beams available in the labora- 



tory, the predictions of the present photon number based 
theory and Glauber photodetection theory are usually 
indistinguishable. 

The number based theory has the advantage that the 
probability is normalizable, for example the probability 
to count one photon in a 1-photon state in the whole 
of space using an array of detectors with rj = 1. The 
Glauber form of the count rate is based on the transi- 
tion probability, however there are advocates for a photon 
number density approach, even within conventional pho- 
ton counting theory. Mandel noted that "there are many 
problems in quantum optics, particularly those concerned 
with photoelectric measurements of the field, which are 
most conveniently treated with the help of an operator 
representing the number of photons" [331 . Mandel and 
Wolf based their general photon counting theory on a 
photon number operator [40,] . Cook observed that there 
is no universal proportionality constant that relates pho- 
ton flux to , and thus the prevailing theory of pho- 
toelectron counting fails to provide a complete descrip- 



tion of photon transport 



He proposed a modified 



photodetection theory based on photon number. A pho- 
ton density ni"' (r, t) , equal to the probability density to 
count a photon at r at time t, is consistent with Cook's 
arguments and with the rules of quantum mechanics. 



VII. ANGULAR MOMENTUM AND BEAMS 

The physical interpretation of the position eigenvectors 
in [Tsj involving AM was motivated by the recent exper- 
imental and theoretical work on optical vortices. These 
vortices are spiral phase ramps described by fields that go 
as exp {ilz^p) and in experiments appear as annular rings 
around a dark center. It can be seen by inspection of ^ 
that the localized states must have orbital AM, and this 
implies a vortex structure that is affected by the choice 



of X- Taking helicity cr = 1 to give a concrete example, 
we can first take to = in ([S]) to give the spherical polar 

vectors {o + ia^^ /%/2 with total AM 0. At 6* = there 

is spin AM h and the orbital AM is —H, while at = tt 
the spin and orbital AM are ~h and h respectively. If 
instead we choose m = 1, the 9 = orbital AM is 0, but 
at = vr it is 2h. For a localized state the vortex has 
not been eliminated, it has just been moved. Thus an 
understanding of optical AM is essential to the physical 
picture of the localized basis states that are used here to 
obtain the photon wave function. 

Theoretically, the simplest beams with orbital AM are 
the nondiffracting Bessel beams (BBs), and these beams 
are closely related to our localized states. They satisfy 
MEs and have definite frequency, cfco, and a definite wave 
vector, kz, along the propagation direction. It then fol- 
lows that the k-space transverse wave vector magnitude 
k± = V^Q ~ ^z' angle 6 = tan^-'^ {k±/kz) also 

have definite values for BBs. Cylindrical symmetry is 
achieved by weighting all (j) equally with a phase fac- 
tor exp {im(t>) . When Fourier transformed to r-space the 
modes go as exp (—ikoct + ilz^p + ikzz) Ji^ {k±r) where 
Iz = m and m ± 1 in ([9]) , J;^ are Bessel functions, ip the 
real space azimuthal angle, and r is the perpendicular 
distance from the beam axis If we select x = so 
that the k-space unit vectors are and 9 in the linear 
polarization basis, B is transverse to z for the 6 mode 
and E is transverse for the (p mode and the linearly po- 
larized modes can be called transverse magnetic (TM) 
and transverse electric (TE) respectively. 

The Bessel functions have a sinusoidal dependence on 
k±r, and this implies that the BBs are standing waves 
that are a sum of incoming and outgoing waves. If inte- 
grated over k± the resulting wave is localized on the z- 
axis at some instant in time that can be defined as i = 0. 
Localization of beams in this way is discussed in [3, . 
If the BBs are then integrated over kz , the result is equiv- 
alent to a sum over all wave vectors, and states localized 
in three dimensions are obtained. But note that this kz 
sum includes waves travelling in the positive and nega- 
tive z-directions. According to the Paley-Weiner theo- 
rem, dkz does not allow exact localization, but this 
restriction does not apply to an integral over all positive 
and negative values. Position is not a constant of the 
motion, and localized states can exist only for an instant 
in time. Exactly localized states in free space are not 
physically possible because they require infinite energy. 
However, our primary concern here is with the use of 
localized basis states for calculation of the photon wave 
function, and we do not require that these basis states 
have a physical realization. 

The real space mathematical description of beams used 
to interpret the AM experiments is usually based on the 
classical energy, linear momentum, and angular momen- 
tum densities. Here, with a basis of position eigenvectors 
in hand that leads to a wave function for a photon in an 
arbitrary state, we are in a position to consider the real 
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space description of the AM of beams from a quantum 
mechanical perspective. The a = 1/2 wave function is a 
solution to MEs, and any derivation based on MEs can 
be adapted to the 1-photon case. The expansion of vec- 
tor potential in [3] that leads to paraxial fields to a first 
approximation can be applied to allow application of our 
formalism to the paraxial beams that are used in most 
optical experiments. Localized states do not exist within 
the paraxial approximation, and the paraxial approxima- 
tion cannot be applied to the position eigenvectors. 

A paraxial beam propagating in the z-direction with 
frequency lo, helicity a, and z-component of orbital AM 
hlz can be described in cylindrical polar coordinates by 
the vector potential [i^ 

A'"*"^ (r, i) = i (x -I- lay) u (r) exp [il^'-p + ikz {z — ct)] . 

(66) 

This vector potential is equivalent to the wave function 
\I>^, (r,t) = (So-,<j'A/2eo/?iA(+) (r,t). The z-component of 
the time average of the classical AM density, equal to 
irx (D*xB + D X B*) , is then found to be 



Jz{r) 



, I . x|2 1 d\u' (r) 
ujh\u{r)\ - 2^^^^ — 



(67) 



separable, the choice x = —4' in Eq. ([7]) gives unit 
vectors (x -I- iay) / \f2 in the paraxial limit which im- 
plies spin quantum number Sz = ±1. The wave func- 
tion (|66p is an eigenvector of Sz with eigenvalue Sz and 
of Lz = —ihd/dip with eigenvalue hlz where f the real 
space azimuthal angle. The latter orbital AM is equiv- 
alent to linear momentum Hlz/r. For this definite he- 
licity state only one term in the photon density (j56p 
contributes. The probability to detect this photon is 
rt'-^-' (r, t) = rS^/"^^ (r, t) , where these field-potential and 
LP densities are essentially equal for a paraxial beam. 
For a coherent state the expansion coefficients Ck.u in the 
1-photon wave function p3p are replaced with the am- 
plitudes 7i(. ^ . Small absorbing particles placed in these 
beams are essentially photodetectors that conserve AM 
by spinning about their centers of mass and rotating 
around the beam axis while they absorb photons 0|. 
The photon number density gives the probability to ab- 
sorb a photon which carries with it spin AM hsz and 
orbital AM hlz ■ For transparent particles the situation is 
more complicated, since re-emission should also be con- 
sidered. 



VIII. CONCLUSION 



It equals the z-component of the AM density ([53]) with 
momentum density given (j51|) or (l52|) without the need 
for time averaging. Thus (p7)) can be interpreted as a 
quantum mechanical AM density that is valid for coher- 
ent and 1-photon states, while ([50|l interpolates between 
these two cases. 

The first term of ([67|) is consistent with orbital AM hlz 
per photon since the photon density given by ([39]) reduces 
ton(i/2) (j.^^) ^ eow|u(r)|Vfi- The last term of ([621) does 
not look like photon spin density. The most paradoxical 
case is a plane wave, as discussed in [i^l- For example a 
wave function proportional to (x -I- icry) exp {ikz — icut) 
implies linear momentum hkz per photon and hence no 
z-component of AM. But we know that such a beam de- 
scribes a stream of photons each with spin AM ha. It was 
observed in 1936 by Beth J8] that a circularly polarized 
beam can cause a disk to rotate, so the beam really does 
carry AM that it can transfer to the disk. The AM of 
this beam resides in its edges, as can be seen from Eq. 
([67| . A new edge is created if the disk intercepts part 
of the beam and this reduces the AM of the beam, al- 
lowing the conservation of total AM [s^ . This is analo- 
gous to the continuum description of a dielectric where 
it is know that the medium is composed of atoms, but a 
continuum description of a uniformly polarized dielectric 
results only is a surface charge. An even closer anal- 
ogy exists between spin AM and a continuous magnetic 
medium where a current in individual molecules reduces 
to a macroscopic current at the edges of the medium. 

In quantum mechanics operators describe observables 
and their eigenvalues are the possible results of a mea- 
surement. While spin and orbital AM are in general not 



We have derived one and two photon wave functions 
from QED by projecting the state vector onto the eigen- 
vectors of a photon position operator. Largely because 
it is still widely believed that there is no position opera- 
tor, this is the first time that a photon wave function has 
been obtained in this way. The two photon wave func- 
tion is symmetric, in agreement with [l| and 0. While 
only the LP wave function gives a positive definite photon 
density, field-like wave functions are widely used and are 
more convenient in many applications. Also, they given 
energy momentum and angular momentum density as in 
([5T|) for example. In the field-like helicity basis the wave 
function pair is 



*(-i/2)(r,i) = 



/fA^'(M), 

-D(+) (r,t). 



(68) 



The wave function components '4'^"^ are given by Eq. 
([33)1 . For definite helicity fields in free space, 6*^+^ and 
the Reimann-Silberstein field vector are just proportional 
to D(+), and thus are equivalent to it. The linear polar- 
ization basis of TM and TE fields can be obtained by tak- 
ing the sum and difference of the definite helicity modes 
as in ([29|) . The photon density is ([56|) 



n, 



(a) 



(r,t) 



c.c. 



(69) 



where n^a^^^ is essentially equal to nS"' except for very 
broad band signals. The 1-photon density can be gener- 
alized to describe the photon density in an arbitrary pure 
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state using the expectation value of the number operator, 

Systematic investigation of photon position operators 
and their eigenvectors clarifies the role of the photon wave 
function in classical and quantum optics. The LP wave 
function defines a positive definite photon number den- 
sity and results in photon wave mechanics equivalent to 
Inagaki's single photon wave mechanics 29]. It is related 
to field based wave functions through a similarity trans- 
formation that preserves eigenvalues and scalar products. 
In free space the field is proportional to the RS 
wave function investigated in [J [3, [l^ |20|. The field 
D(+' (r, t) is proportional to the Glauber wave function 
[1, H, which gives the photodetection amplitude for 
a detector that responds to the electric field ll|. While 
only fields and potentials are locally related to charge and 
current sources, Fourier transformation of k-space prob- 
ability amplitudes naturally leads to the LP form [a, [13] . 
The similarity transformation between the field-potential 
and LP wave functions makes the choice a matter of con- 
venience for most purposes. 

By the general rules of quantum mechanics the LP 
wave function is the probability amplitude to detect a 
photon at a point in space. It and the closely related 



field-potential wave function pair obtained by solution of 
MEs are ideally suited to the interpretation of photon 
counting experiments using a detector that is small in 
comparison with the spatial variations of photon density. 
It is not subject to limitations based on nonlocalizability, 
and coarse graining or restriction to length scales smaller 
than a wave length is not required. Exact localization 
in vacuum requires infinite energy and is not physically 
possible, but position eigenvectors provide a useful math- 
ematical description of photon density . Photon number 
density is equivalent to integration over undetected pho- 
tons in a multiphoton beam. In an experiment where ab- 
sorbing particles are placed in a beam, the particles act 
as photodetectors which can sense the spin and orbital 
angular momentum of the photons. Our formalism justi- 
fies the use of positive frequency Laguerre-Gaussian fields 
as photon wave functions and gives a rigorous theoreti- 
cal basis for extrapolation of their range of applicability 
from the many photon to the 1-photon regime. 
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